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METHOD  OP  CALCULATION  OF  NONSTATIONARY 
AERODYNAMIC  CHARACTERISTICS  OF  A LOW-FLYING 
WING  WITH  A CYLINDRICAL  FUSELAGE 


L.G.  Tsvetkov 

Department  of  Applied  and  Computer  Mathematics 

The  article  gives  an  account  of  the  method  of  calculation 
of  nonstationary  aerodynamic  characteristics  of  a low-flying 
wing  of  complex  configuration  with  a fuselage  In  the  form  of  an 
infinite  cylinder  taking  into  account  the  nonlinearity  caused  by 
the  effect  of  the  setting  angle  of  attack  of  the  wing  and  the 
movement  of  points  of  the  surface  of  the  wing  with  respect  to  the 
wing  with  nonstationary  motion  of  the  complex. 

Let  us  examine  the  motion  of  a wing  with  a fuselage  in  a 
nonviscous  incompressible  fluid  with  a constant  speed  Uq  in 
parallel  to  the  reference  surface.  Let  us  assume  that  the  com- 
plex accomplishes  harmonic  vibrations  with  low  relative  freauency 
made  up  from  the  translational  vertical  movement  and  rotation 
around  the  OX  and  OZ  axes  of  the  bound  coordinate  system  OXYZ 
(see  the  figure).  Vibrations  of  the  complex  are  accomplished 
with  respect  to  the  average  position  at  which  the  axis  of  the 

fuselage  is  parallel  to  the  reference  surface  and  is  distant  from 

• • 

it  by  the  magnitude  , where  B Is  the  root  chord  of  the 

wing. 

Let  us  represent  the  thin,  slightly  curved  wing  of  complex 
configuration  in  the  form  of  a combination  of  plates  and  plans, 
and  let  us  assign  the  geometry  of  the  complex  by  means  of  a,  - 
aspect  ratio,  ij,-  taper,  setting  angle  of  attack, 


sweep  angle  of  the  leading  edge,  Y* - deadrise  angle  of  the  kth 
nlan  of  the  wing,  e="J"~  rel&tive  rise  of  the  point  of  inter- 
section of  the  root  chord  above  the  axis  of  the  fuselage  and 

R = -| — the  relative  radius  of  the  fuselage. 


Let  us  present  the  coefficients  of  aerodynamic  force  and 
moments  of  the  lifting  wing  by  coefficients  of  rotating  deriva- 
tives in  the  form 


S = V2I  \ V , m,=iV‘£mlv  , rn.rm’o  ♦ , 

v=<p,h,o,<p,fi.c  , 


(I) 


where  is  the  small  increase  in  the  trim  angle  of  the 

complex,  >i=fc(t}-  increase  in  the  flight  altitude,  ■c-  = w(t) 
small  increase  in  the  bank  angle. 

The  streamline  flow  around  the  fuselage,  which  has  the  form 
of  an  infinite  cylinder,  vibrating  above  the  screen,  is  modeled 
by  the  system  of  dipoles  with  a constant  (along  the  OX  axis)  in- 
tensity with  vertical  vibrations  and  with  an  intensity  changing 
according  to  the  linear  law  with  rotating  vibrations  with  resnect 
to  the  OZ  axis.  The  position  of  the  dipoles  and  their  intensity 
are  determined  by  the  principle  of  inversion  with  respect  to  the 
cylindrical  surface  and  the  principle  of  mirror  image  with  respect 
to  the  screen,  assuming  in  the  first  approximation  the  intensities 
of  the  dipoles  to  be  the  same  as  those  with  vibrations  of  the 
cylinder  in  an  infinite  medium, [1]. 


We  find  the  speeds  caused  by  vibrations  of  the  inherent 
fuselage  above  the  screen  at  point  • Cj  of  the  bound 

coordinate  system  by  the  method  of  successive  approximations 
according  to  the  following  equations  in  which  all  the  linear 
dimensions  belong  to  the  root  chord  of  the  wing: 
u,(t)  - with  translational  vertical  vibrations  with  velocity 

w*r°*  <s) 
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with  rotational  vibrations  with  the  angular  velocity 

. <3> 


where 


'■r%£  tHv.*  • 

i 

J>i=Vv -)*♦*!  • 


The  quantities  m»  and  *}„  entering  into  the  velocity 
expressions  are  connected  by  the  dependences 


_T_  T 


T*i  mT*4 


* — nii  » t /«.  t,  r»*  u 

* ’ « ’TtT?  ’ 7l*  = -(zh.4Tl«)>  7«  =^T  • 


In  the  first  approximation  (y  = 1)  we  must  assume  that 

m*«  I.  n',-  0. 

We  model  the  vortex  surface  of  the  wing  by  the  system  of 
oblique  horseshoe-shaped  vortices,  and  we  assume  that  the  free 
vortices  are  are  first  located  in  the  plane  of  the  wing  and 
behind  the  trailing  edge  according  to  the  velocity  of  the  non- 
disturbing advancing  flow  in  parallel  to  the  reference  surface. 
According  to  recommendations  of  work  [2],  we  produce  the  location 
of  the  horseshoe-shaped  vortices  on  each  plan  of  the  wing  and  the 
selection  of  the  check  points  in  which  the  boundary  conditions  on 
the  wing  are  satisfied.  The  numbering  of  the  check  points  is 
continuous  for  the  whole  right  half  of  the  lift  wing;  the  reading 
is  conducted  in  each  plan,  beginning  from  the  first,  along  the 
span  from  the  root  to  the  end  chord  of  the  plan  and  along  the 
chord  from  the  leading  edge  to  the  trailing  edge,  where  is 

the  number  of  vortices  located  along  the  span  of  the  plan,  n^  - 
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the  number  of  vortices  along  the  chord. 


To  fulfil  the  boundary  condition  of  sealing  on  the  reference 
surface,  the  - vortex  system  of  the  wing  is  mirror  reflected  with 
resnect  to  the  plane  of  the  screen  with  opposite  signs  of  the 
vortices  . 

The  fufllment  of  the  condition  of  sealing  on  the  cylindrical 
fuselage  with  retention  of  the  boundary  condition  on  the  refer- 
ence surface  is  ensured  by  the  successive  inversion  of  each  1th 
transformed  horseshoe-shaped  vortex  and  its  image  relative  to  the 
screen  with  respect  to  the  cylindrical  surface  and  mirror  image 
relative  to  the  reference  surface  [3].  The  transformed  horseshoe- 
shaped vortex  is  obtained  by  the  projection  of  the  initial  i th 
vortex  of  the  wing  onto  the  plane  which  is  parallel  to  the  axis 
of  the  fuselage  and  passes  through  the  trailing  edge  of  the  plan. 
As  a result  of  a similar  image  within  the  fuselage,  there  is 
available  a number  of  horseshoe-shaned  vortices  of  the  same  in- 
tensity as  that  of  the  ith  vortex  of  the  wing  but  having  dif- 
ferent parameters.  Replacement  of  the  initial  horseshoe-shaoed 
vortex  by  the  transformed  vortex  with  use  of  the  method  of  suc- 
cessive inversion  is  admissible  in  view  of  the  small  setting 
angles  of  attack  of  the  wing. 

Parameters  of  the  inherent  ith  transformed  vortex  ( - 

i 

semispan,  - deadrise  angle,  \i  - sweep  angle  of  the  connected 

vortex,  • ’ll  £;  - coordinates  of  the  position  of  the  vor- 

tex in  the  bound  coordinate  system)  are  calculated  in  terms  of 
the  initial  geometric  parameters  of  the  complex  in  the  following 

“ay:  , m«m.) 


where 


* ".n„,  1 

» _ 

A 

2-n  ii 

p»t  .i; 

X t *(<• 

ms 

| VR‘-e1' 

| e | » R 


y (4) 


when  k = 0 we  must  assume  that  ■ 0 <u~J-  X,  operation  of 

the  separation  of  the  fraction  part;  ()  - operation  of  the 
separation  of  the  whole  part  of  the  number. 


The  conformity  between  the  number  of  the  vortex  i and  the 
number  of  the  panel  k,  on  which  this  vortex  is  located,  is  found 
from  the  condition  of  the  fulfilment  of  the  equality 


T~  u N <«i  « Hn.M,  . 

P»«  r r P'0 


For  an  account  of  the  change  in  velocities  at  the  check 
points  j of  the  surface  of  the  wing  owing  to  movements  of 
the  latter  with  respect  to  the  vortices  reflected  relative  to 
the  screen  with  nonstationary  motion  of  the  complex,  the  veloci- 
ties induced  by  these  vortices  we  expand  in  Taylor  series  with 
respect  to  the  small  parameters  of  the  movements  % and 
and  let  us  be  limited  to  terms  of  the  first  order  of  smallness, 
so  that  in  the  case  of  the  constant  (in  time)  circulation  of 
the  horseshoe-shaped  vortex  ( r^lt)  -u,B  r -const  ) the  value 

of  the  velocity  is  determined  by  the  equation 


AC; 


i »C 


(5) 


With  the  harmonic  law  of  the  change  in  the  circulation  of 
the  vortex  with  time,  the  increase  in  velocities  at  points  of 
the  surface  of  the  wing  owing  to  their  small  movements  has  a 


higher  order  of  smallness,  and  the  value  of  velocities  is  deter- 
mined by  a well-known  manner  (at  values  of  the  Strouhal  numbers 

PS 

tending  to  zero,  ^ = — 0); 

- with  the  sinusoidal  law  of  the  change  in  circulation  of 
the  vortex  (I^.su^SItsun  pt) 


- with  a change  in  the  circulation  according  to  the  cosine 

law  yt) 


pt*  Vos 


With  the  harmonic  law  of  the  change  in  kinematic  parameters 
of  motion  of  the  complex 

f (t)  = f stn  P,  t , S (t)  = K sun  pjt  t ©(t j = ©* sun  p^t 

the  increase  of  the  coordinates  of  the  check  point  In  the  coordi- 
nate system  connected  with  the  horseshoe-shaoed  vor?extref¥ec?ed 

* 

relative  to  the  reference  surface  have  the  form 

♦ 2h*cos¥i>«-Tipzt~2o*[^tosyi»(iji«h,)sinyi1sinpil^costti  , (8) 

2ft* sun  sin  p2t  - 

“ 2»*[CjSun  ♦ flouts  ipjsun  p4t . 

the  increase  in  the  coordinates  of  the  cheek  point  in  the 
coordinate  system  of  the  vortices,  which  are  a mirror  image 
with  respect  to  the  screen  of  the  vortices  located  within  the 
fuselage,  has  a simple  form 

Alls(f *4 sln Pi1  ♦ P,t)cosy J-  o*(5iCosiji^ , , 

A>IJ= *n  PA  * 2 *’  ptt)sin  Y l - sun  ij . tos  ) sun  Pjt . 
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The  coordinates  of  the  check  Doints  entering  into  these 
expressions  are  comouted  in  terms  of  the  initial  geometric  param. 
eters  of  the  complex 


* Ip 

ft  - f 

V = ' 6 + £ l*  stn  Vp  + /*i f i - i ^ -^) tw  fj ^ f 

r-4 

cos  if.  +(lK-l-)  sin  f . ■ oci 


In  turn 


where  and  jj  are  computed  according  to  equations  (4) 

with  the  replacement  of  i by  y . 


Representing  the  dimensionless  circulations  of  the  vortices 
in  the  form  similar  to  the  expansion  of  (1) , using  the  dependences 
(2),  ( 3 ) , (5),  (6),  (7),  (8),  and  (9)  and  satisfying  at  the 
check  points  the  boundary  condition,  which,  taking  into  account 
the  nonlinearity  according  to  the  setting  angle  of  attack,  has 
the  form 

W^-u.sLnocj  cosoc-cds^sLti p^t  tptU* cosoc^os^j  cosp^  + 

+ P«5'AV’nV  SitoseN  C0S  1,i)C0&P<t*  - sLn  f 

we  obtain  a number  of  systems  of  linear  algebraic  equations  for 
determining  the  intensities  of  the  horseshoe-shaped  vortices  of 
the  wing 
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fc«l  « 

£r>  =-2Cr*tVVV’  <•  I 

k*l  I k»4  4 4 

gr’Fs.i^rtl(^.»4*M-S^*(li*‘A‘lH)‘Wl»%il ' 

E r’ F;  = **  (ijtj* j * 5j  “Mr  li ) ■ - s r‘  s - 

la)  4 

i=  4 4 

£ r* hi r - V3r CBS  fi 4 ^ 5111  f’i ) ' S r ^ a li  * 

i * 4 « 

j =<.2,3 m = £>A., 


where 


F..=f..+irf.  -f..-H.. , 

kt  Jii  Jli 

=i..CDS«j.i^i..COS%  , 
&i.=pt.CQS,Pi*&pi.cOS  ?<}>;, 

VV^fV^i  » 

B^o’.cos^-e;.^^;. 

VV'S. 


VV'VV^ ; 

&ii=Piislnfi+5'^islTl?9i.  > 

vv^’V^v 

Lis*;i“5ti“»,“1sinifi  ■» 


S =V.  ■►  frS: . 
ll  ll  U 


In  turn 


L 


* 

, t0S(0t;-OC.) 

fcf  ^(fi  - *fl)  “tosi TT-  - (t . - Yl)  - 

" C [»^os(«Pj  -flv  " 4>[ )] + |l  uV>tn(fJ-  Tj], 

flu},  , nos(oc;-«.)  Bltf, 

-Zl[^  tns(ti  - )-'^ )]♦  c05Wr^)-^  , 

b^=  £ [(6-^  «»»*[♦  ^ *«.  )c°Hf j - ^)  ■ (^ CDS  t I + T;  Sin^  )sln(v1  - *f  f)]-  1 

- -T  "Bixs-Z  _i\  / -r.  <ih)J  ~r  -i\  • , -T\l 

V C [( 0^ stn r i - af  cosf  JC0SWj- sinVL - if “** Pin  Wr? i >J ' 

- C [( g^j*  S1-n  fi  - ^ ™ f i )cosCvi  '4>i ) - ( 0^  stn H>c  - -jj-  COS  Vi  )unty  - ^ ) j , 

2i»a  , - .COS^.-K.^  flu  , _ , 

“(trfil-s*-  * jf“"(v*). 

We  determine  the  quantities  . jk.  .rijL.,  ^k.  etc.  by  the 

same  equations  (11),  having  replaced  the  parameters  of  the  i th 
vortex  |k.  •1Kj  *Sij  • Vi  • \l  by  parameters  of  the  vortices 

and  so  on. 


The  parameters  for  the  calculation  of  velocities  induced 
by  the  intrinsic  i th  vortex  of  the  right  half  of  the  wing,  the 
vortex  symmetric  to  it  relative  to  the  wing,  and  also  vortices 
of  the  left  half  of  the  wing  can  be  determined  by  coordinates  of 
the  check  Doint  and  parameters  of  the  transformed  vortex  in  the 
following-  manner: 


■■iir 


Ti.i=t1?r1/^C0Sv,i-^r^^slTl,i’i3CDSOCi4*» 

7i/[(l?J<',/i+2flo)C0Syi+^i-tl)^n'('i]t05<il- X,  iif-tyi+H i^sta^Kj-CJrasfi , 
^‘^'^“•TiHCj+CiJsinvJcoscii**,  t^)tDiVi , 

ff’iif[^i*7i+2,k«)C05Vi-(Cj*5t)»l,lTjo>«r*,  fffii=('lit1i*2R.)sinfiHCi^i)c'>s'Pi: 

*=(V  Vsit,oci  • 


Parameters  for  the  calculation  of  velocities  under  the  sum 
sign  with  resoect  to  y are  determined  by  equations 


Fr  -c1  — B"F^ r H £ T . aT  = \y=-yxT.=n7C  si.n  ■■  ■ » 

S S N S'  WV\i*W 

«^=-<j>T=-lfyj^&(j}7=s-ausin  ~^fi  = iV^W^T  » 

^=(Va¥)“s^'(Ci"t"^)sLl’^ » (ci‘  ^)tK^fc'qT'^in^' 

^®(Cj “^r“)n,5ti“ (if  ^T" )*n  ^ » 

rVfc'^T')ws^(ci+  ^=fej+^W‘lv  nH611^ » 


>(I2) 


where 


=T_x 


i to  X 


ftox 


r*<  5 R S R r* i t4R 

” = =(ty,(dV  V)^  rj’ 


dr  - _ji  T*  1 


(fi*)MdT)x 

ir=-(2Ka+Qr),  cT=c\  8T=-U»«.+  8r),  dr=dr. 


(13) 


We  compute  the  quantities  a',  b',  c',  and  d’,  which  corres- 
Dond  to  the  first  approximation  (Y  = 1),  according  to  equations 
(13),  having  assumed  here  that 

We  compute  the  parameters  for  the  calculation  of  velocities 
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under  the  sum  sign  with  respect  to  X according  to  the  same 
equations  (12),  having  assumed  in  the  first  approximation  (X 


that : 


We  compute  the  velocities  ,avw*"  an<^  derivatives  of 

velocities  according  to  theStrouhal  number  according  to  equations 

of  work  [4],  producing  preliminarily  a demeasuring  by  means  of 

the  root  chord  of  the  wing;  derivatives  of  velocities  it*- . . ^r3- 

~ i . H../1S  •’  *c  *> 

can  be  obtained  by  means  of  differentiating  with  respect  to  n and 

£ of  expressions  for  the  appropriate  velocities. 


After  solving  the  system  of  algebraic  equations  (10),  it  is 
possible,  by  using  the  Zhukovskiy  formulas,  to  determine  the  to- 
tal aerodynamic  characteristics  of  the  lifting  wing  with  a 
cylindrical  fuselage: 

where  S is  the  area  of  the  projection  of  the  cantilevers  onto 
the  plane  XOZ  determined  by  equation 


ii 


Bibliography 


1.  Kudryavtseva,  N.A.  Streamline  flow  of  an  infinite  cylinder 
carrying  out  vibrations  in  an  incompressible  fluid.  Trudy  TsGI, 
Issue  810,  1961. 

2.  Belotserkovskiy , S.M.  Thin  lifting  surface  in  a subsonic 
gas  flow.  Nauka  , 1965. 

3.  Tsvetkov,  L.G.  ADoroximation  method  of  the  calculation  of  the 
distribution  of  the  aerodynamic  load  along  a low-flying  wing 
with  a fuselage.  Trudy  LKI  , Issue  XIX,  1070. 

4.  Kolesnikov,  G.A.  Three-dimensional  velocity  field  of  non- 
stationary oblique  vortex.  Trudy  TsAGI,  Issue  810,  1961. 


13 


DISTRIBUTION  LIST 


Bfc 


DISTRIBUTION  DIRECT  TO  RECIPIENT 


ORGANIZATION 

MICROFICHE 

ORGANIZATION 

MICROFICHE 

A205 

DMATC 

1 

E053 

AF/INAKA 

1 

A210 

DMAAC 

2 

E017 

AF/  RDXTR-W 

1 

B3UU 

DIA/RDS-3C 

8 

E404 

AEDC 

1 

CO  1*3 

USAMIIA 

1 

E408 

AFWL 

1 

C509 

BALLISTIC  RES  LABS 

1 

E410 

ADTC 

1 

C510 

AIR  MOBILITY  RAD 

1 

E413 

ESD 

2 

LAB/FIO 

FTD 

C513 

PICATINNY  ARSENAL 

1 

CCN 

1 

C535 

AVIATION  SYS  COMD 

1 

ETID 

3 

•gppi 

■■H 

t 

NIA/PHS 

1 

C591 

FSTC 

5 

NICD 

5 

C619 

MIA  REDSTONE 

1 

D008 

NISC 

I 

H300 

USAICE  (USAREUR) 

1 

P005 

ERDA 

X 

P055 

Cl A/CRS /ADD/SD 

1 

NAVORD8TA  (50L)  1 . 

NASA/KSI  1 

AFIT/LD  1 


